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Brief History of Related Events 



Historical Perspective 

February 2000 
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Absolute Radiometric Calibration 


The purpose of absolute radiometric calibration is to 
determine the expected response for any in-band entrance 
aperture radiance presented to the telescope. 


CalCoef : 



where: CalCoef = In-Band Calibration Coefficient (DC/(W/m 2 -sr)) 
DCj = Measured Digital Counts in the j th band (DC) 

Lj = Entrance Aperture Radiance in the j th band (W/m 2 -sr) 
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Overview of Original Method Used to Establish 
Absolute Radiometric Accuracy of Remote 
Sensing Instruments Using Stellar Sources 



Original Method to Establish Absolute Radiometric 
Calibration Accuracy Using Stellar Sources 


■ Identify Acceptable Radiometrically Characterized Stars 

— Compute Broadband Entrance Aperture Flux Density 

■ Gather Stellar Response Data 

— Extract Stellar Scene Data (sum digital counts) 

■ Generate Calibration Coefficient 

■ Verify Results Using Independent Sources 

■ Determine System Stability 

■ Determine System Linearity 
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Bright Star Catalogue Stars 
-1 .46 < Vmag < 8.0 = 91 1 0 Items 
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Bright Stars most likely to be compatible with 
IKONOS-like Remote Sensing Instruments 


Bright Star Catalogue Revision 5.0 Histogram 
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BSC Remote Sensing Compatible Stars ( ) 

and Spectrophotometric Standard Stars (28 
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Remote Sensing Compatible Spectrophotometric 

Standard Stars (28 (Excluding Vega and 1 repeated star) 


RS ID 

BSC ID 

RA (J2000) 

DEC (J2000) 

Proper Name 

Vmag 

Spectral Type 

HST 1 

HR7001 

18:36:56 

38.7836 

Alpha Lyrae 

-0.03 

A0V 

HST 2 

HR718 

2:28:09 

8.4601 

Xi 2 Ceti 

4.28 

B9III 

HST 3 

HR1544 

4:50:36 

8.9002 

Pi 2 Orionis 

4.36 

A1V 

HST 4 

HR3454 

8:43:13 

3.3986 

Eta Hydrae 

4.30 

B3V 

HST 5 

HR4468 

11:36:40 

-9.8023 

Theta Crater 

4.70 

B9.5V 

HST 6 

HR4963 

13:09:56 

-5.5390 

Theta Virginis 

4.38 

A1IV 

HST 7 

HR5501 

14:45:30 

0.7174 

108 Virginis 

5.68 

B9.5V 

HST 8 

HR7596 

19:54:44 

0.2735 

58 Aguilae 

5.62 

AON 1 

HST 9 

HR7950 

20:47:40 

-9.4958 

Eta Aguarii 

3.78 

A1V 

HST 10 

HR8634 

22:41:27 

10.8314 

Zeta Pegasi 

3.40 

B8V 

HST 11 

HR9087 

0:01:49 

-3.0275 

29 Piscium 

5.12 

B7III 

HST 12 

HR153 

0:36:58 

53.8969 

Zeta Cassiopeiae 

3.66 

B2IV 

HST 13 

HR1996 

5:45:59 

-32.3065 

Mu Columbae 

5.17 

09V 

HST 14 

HR4554 

11:53:49 

53.6948 

Gamma Ursae Majoris 

2.44 

A0V 

HST 15 

HR5191 

13:47:32 

49.3133 

Eta Ursae Majoris 

1.86 

B3V 

G&S 1 

HR6787 

18:08:45 

20.8144 

102 Her 

4.36 

B2IV 

G&S 2 

HR3454 

8:43:13 

3.3986 

Eta Hydrae 

4.30 

B3V 

G&S 3 

HR6588 

17:39:27 

46.0064 

85 lot Her 

3.80 

B3IV 

G&S 4 

HR6556 

17:34:56 

12.5600 

55 Alp Oph 

2.08 

A5III 

G&S 5 

HR6561 

17:37:35 

-15.3986 

55 Xi Ser 

3.54 

FOlVDel Set 

G&S 6 

HR1411 

4:28:34 

15.9622 

77 Thel Tau 

3.84 

K0II lbFe-0.5 

G&S 7 

HR4247 

10:53:18 

34.2150 

46 LMI 

3.83 

K0+III-IV 

G&S 8 

HR4335 

11:09:39 

44.4986 

52 Psi UMa 

3.01 

Kim 

G&S 9 

HR5854 

15:44:16 

6.4256 

24 Alp Ser 

2.65 

K21 1 IbCN 1 

G&S 10 

HR3905 

9:52:45 

26.0069 

24 Mu Leo 

3.88 

K2IIICN1Ca1 

G&S 11 

HR3845 

9:39:51 

-1.1428 

35 lot Hya 

3.91 

K2.5 1 1 l-l 1 1 bBa0.3 

G&S 12 

HR2574 

6:54:11 

-12.0386 

14 The CMa 

4.07 

K4III 

G&S 13 

HR3003 

7:46:07 

18.5100 

81 Gem 

4.88 

K4 1 1 1 - 1 1 1 b 

G&S 14 

HR2905 

7:35:55 

26.8958 

69 Ups Gem 

4.06 

MOIII-lllb 

G&S 15 

HR1556 

4:52:32 

14.2506 
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Remote Sensing Compatible Spectrophotometric 
Standard Stars Spectral Data 
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Hertzprung-Russell Diagram Showing Remote 
Sensing Characterized Star Categories 
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Dec (deg) (J2000) 


Brightest G&S Stars and Those Used for IKONOS 
Absolute Radiometric Accuracy Assessment (7) 
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Flux (mW/m z -Hz) 


J.E. Gunn & L.L. Stryker - Stellar Spectrophotometric Atlas 
3130 < X < 10800 A (1983), Remote Sensing Stellar Subset 



Frequency (Hz) 
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Un-normalized Stellar Spectral Data v. Frequency 

Remote Sensing Stellar Subset 
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Flux (mW/m 2 -pm) 


Un-normalized Stellar Spectral Data v. Wavelength 

Remote Sensing Stellar Subset 
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Flux (mW/m 2 -|jm) 


De-reddening Coefficients for Stellar Spectral Data 

Remote Sensing Stellar Subset 



Wavelength (pm) 
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Flux (mW/m 2 -pm) 


Re-reddened Stellar Spectral Data v. Wavelength 

Remote Sensing Stellar Subset 



Wavelength (pm) 
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Flux (mW/m 2 -gm) 


Subset of G&S Stellar Spectra Used To Assess 
IKONOS Absolute Radiometric Accuracy 



Wavelength (gm) 
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IKONOS MS Bands Relative Spectral Response (RSR) 

(http://www.spaceimaging.com/products/ikonos/spectral.htm) 


Lower 50% Upper 50% Bandwidth Center 


Band (nm) (nm) (nm) (nm) 


Pan 

525.8 

928.5 

403 

727.1 

MS-1 (Blue) 

444.7 

516.0 

71.3 

480.3 

MS-2 (Green) 

506.4 

595.0 

88.6 

550.7 

MS-3 (Red) 

631.9 

697.7 

65.8 

664.8 

MS-4 (VNIR) 

757.3 

852.7 

95.4 

805.0 


Note: 

1. Spectral Bandwidths are Full-Width at Half-Max (FWHM) 

2. Panchromatic Band does not incorporate spectral filtering - response is that of 
optics/detector only 


Values Shown in Color are Bandpass Definitions Using 50nm-interval Datasets Without Interpolation 
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IKONOS Multi-spectral Bands RSR - 50nm Intervals 
(as published on GeoEye website) 





Flux (mW/m 2 -pm) 


G&S Stellar Spectra and IKONOS MS Bands RSR 



Wavelength (pm) 
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Digital Counts 


IKONOS Stellar Calibration Coefficients 2001 
Compared with JACIE 2001 Results 
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IKONOS Stellar Calibration Coefficients 2001 
Compared with JACIE 2001 Results 
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Considerations to Improve the Stellar 
Calibration Approach 



Traditional Remote Sensing System Calibration 
Methods Determine Entrance Aperture Radiance 
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Typically the System Response Function Takes Into 
Consideration All of the Components of the System 
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Calibration Coefficient is the Relationship between 
Entrance Aperture Radiance and Digital Counts 



Digital Counts ( DC 
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' w \ 
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Stellar Calibration Uses Well-characterized Stellar 
Flux Densities Arriving at the Payload Aperture 


Digital 

Counts 

(DC) 
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Stellar Calibration Method as Presented in 2002 
ISPRS Paper 

■ Assumed 

— All of the energy was represented as if it was being captured within a 
single Instantaneous Field of View (IFOV) of the detector 

■ Did not consider Encircled Energy (s f ) function for the optical system 

— Encircled Energy Function determines the amount of energy falling on the 
focal plane that is higher than the response threshold of a detector circuit 

■ Did not consider compression artifacts imposed on stellar scene data 

— Older Compression Algorithms and Points Sources don’t play well 

■ Did not include an absolute spectral reference associated with the 
vicarious ground calibration approach 

— Stellar Spectral Content versus Solar/2xAtmosphere/Ground Spectral 
Content 
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< < 


Pixels with Measurable Response Summed to 
Energy Contained Within Point Spread Function 




Average Number of Detectors within 7x7 pixel Area 
Summed to Produce the Total Digital Counts 
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Stellar Calibration Sensor Response Function 
Algorithm 



"oj-p/n 

J -^ J -*RSR,Ws, i (A m )*SCF IJ *dX 

350 nm do 


where: Lj pd) = Equivalent Radiance of Star at the Aperture (W/m 2 -sr) 
F|(X) = Spectral Flux Density at Entrance Aperture (W/m 2 -pm) 
Q d0 = Solid Angle of a single Detector (sr) 

RSRj(X) = Relative Spectral Response of jth Band (%) 

Sfj (A m ) = Encircled Energy Function of jth Band (%) 

SCFj j = Spectral Correction Function of jth Band (%) 
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CRSS w/simulated stellar PSF (log scale) 

15x15, 48um pixels, X = 0.75 um, Focus at 700 km 


1 E + 0 


ill 


fEiimassam! 

ee Sf fil e||.3e 

Mil® 

■-0"" 

■ii i 

BH 


3BBBMBBHDffly|piWll 

FFt™=n smt a™ mas sm jam mm (WT* HHH1 H I IH H 1 1 




SHIIIIllllll' 

w * * ™ w §« Hi Sflffl SSSiffi:: 





i E-i 


IE-2 


1 E-3 


1 E - H 


1E-E 


HHIllllHlIlim pfTT 


1 E-6 


LOG HUYGENS PSF 


CRSS , 7 FLOS, REG RPE , MOVE Fill -Z 0-2, RDD 0,5 DEG FM2 
SHT MRR Hi 2006 

,7500 Tn 0,7500 //m RT 0,2330, 0,5000 EEC, 

I URGE SIZE IS 768 , 00 //m SQURRE 
STREHL RRTIO 0,799 


CENTER COOEEHNRTES: - 4 , 076 1 063 IBS 


1 , 939982H6E+00 1 


C R S S *P X M R , ZMX 
CONFIGURRTION 1 OF ll 


ITT 


2006 JACIE Stellar Calibration March 15, 2006 


35 






Point Source Illumination at Focal Plane Is 
Distributed According to Optical Transfer Function 



Detectors 
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Encircled Energy Function Describes Unaccounted 
Energy Due to Spatial Distribution at Focal Plane 
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Discrete Detectors in the Focal Plane Under- 
sample the Point Source Energy Distribution 
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Low-contrast Stellar Scenes Cause Compression 
Artifacts that Contribute to Algorithm Uncertainty 
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Encircled Energy Function 


Empirical Encircled Energy Functions that Provide 
Match Between Stellar and Ground Results 
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Spectral Correction Function can be Employed to 
Reference the Stellar to a Solar/Ground Spectra 


1 100 nm 

\L si jmSRU)ca 

350 nm 


SCF i = 


\L sta MW 

4 

1 lOO^m 

lL gHd U)RSR(A)cU 

350 nm 



where: L star (X) = Equivalent Radiance of Star at the Aperture (W/m 2 -sr) 

L gnd (l) = Radiance of Reference Ground Scene at the Aperture (W/m 2 -sr) 
RSRj(X) = Relative Spectral Response of jth Band (%) 

SCFj j = Spectral Correction Function of jth Band (%) 
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Correlation between Corrected Stellar Calibration 
and Vicarious Ground Calibration Results 
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Topics for Future Investigations 


■ Add New Characterized Stars to Collection Plans 

■ Revisit Stellar Imaging Scenario 

— Develop Encircled Energy Function 

- Linear Array 

- Push-broom Scanning Systems 

— Separate data from compression noise 
— Determine uncertainty of measurement 

■ Collect and Analyze Additional Stellar Imagery 

— Analysis primarily uses only 7 stellar collects 
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Summary 


■ Additional research needed to further refine the algorithm 


* 


ITT 


2006 JACIE Stellar Calibration March 15, 2006 


44 




